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Analytic solutions for the three-gluon and ghost-gluon vertices in Landau gauge Yang-Mills theory at 
low momenta are presented in terms of hypergeometric series . They do not only show the expected 
scaling behavior but also additional kinematic divergences when only one momentum goes to zero. 
These singularities, which have also been proposed previously, induce a strong dependence on the 
kinematics in many dressing functions. The results are generalized to two and three dimensions and 
a range of values for the ghost propagator's infrared exponent k. 
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I. INTRODUCTION 

During the last ten years the infrared (IR) behavior 
of Yang-Mills Green functions was studied intensively. 
Especially the propagators are of interest as they are 
linked to the confinement problem via the scenarios of 
Kugo-Ojima [1] and Gribov-Zwanziger j2]|3]- Both pre- 
dict in Landau gauge an IR enhanced ghost propaga- 
tor that is responsible for the appearance of long-range 
forces, whereas the gluon propagator should vanish (or 
be at least finite for the former). Indeed such behav- 
ior was found by functional methods (Dyson-Schwinger 
equations (DSEs) g] |5] EJ El E2 as well as Renormal- 
ization Group (RG) equations |9l ITO]). 

The DSEs of Yang-Mills theory in Landau gauge allow 
two different kind of solutions, called the scaling and the 
decoupling solution in ref. [9]. The latter has been ad- 
dressed in refs. [9] EU H] [13] [H E] ■ Especially, for this 
solution the behavior of vertices was determined in [T3] 
as trivial, i.e. the dressing functions are at most constant 
in the IR, and there is no IR enhancement. Further- 
more, the decoupling solution does break global BRST 
and global gauge invariance in the confining phase j9], 
whereas the scaling solution can provide a picture for 
confinement in Landau gauge |16l IT?] . Therefore we will 
concentrate on the scaling solution, as described in sec. 
II, in this paper. 

Besides the propagators also the vertices play an 
important role in Landau gauge Yang-Mills theory. 
Whereas the ghost-gluon vertex is IR regular the three- 
gluon vertex is strongly IR divergent [18]. In fact, this 
happens for all dressed gluonic vertices in the uniform 
limit that all momenta vanish simultaneously. Thereby 
the IR behavior of any diagram in a DSE is characterized 
by its bare vertex [19] . A recent power counting analysis 
showed that the three-gluon vertex can feature additional 
kinematic singularities when only a single gluon momen- 
tum vanishes [13]. As has been explained in detail in 



ref. [T7], the importance of the vertex functions is even 
more pronounced in the quark sector: Here the dressed 
quark-gluon vertex may provide an explicit mechanism 
for quark confinement in Landau gauge QCD. In this ver- 
tex it is a strong soft-gluon divergence in quenched QCD 
that is able to generate the necessary interaction strength 
to permanently confine quarks. This strong divergence 
is in turn self-consistently induced by a contribution in- 
volving effectively the strongly IR divergent three-gluon 
vertex. 



Here we present a detailed study of the IR behav- 
ior of the three-point vertices (the three-gluon and the 
ghost-gluon vertex) in Landau gauge Yang-Mills theory. 
Within a semi-perturbative approximation we provide 
explicit analytic results for the full kinematic dependence 
in the IR limit that confirm and exemplify the results of 
our previous power counting analysis [Tjt]. These results 
allow to develop better ansatze for these vertices in sub- 
sequent investigations. This will provide more reliable 
results for the gluon propagator and the full momentum 
dependence of Yang-Mills vertices and in particular for 
the analysis of the quark-gluon vertex. 
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We first review the IR fixed point structure of Lan- 
dau gauge Yang-Mills theory obtained from a pure power 
counting analysis in sec. [H] To go beyond this we present 
an analytic approximation scheme that covers the lead- 
ing IR contributions to the three-point vertices in sec. 
|III| where we also present general results for the trian- 
gle graphs with dressed propagators in the IR limit and 
show how they can be evaluated in the Euclidean case. 
The explicit results for the individual tensor structures 
of the three-gluon vertex and the ghost-gluon vertex are 
presented in sees. IV and [V] respectively. Technical as- 
pects like the tensor decomposition are deferred to two 
appendices. 
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II. INFRARED BEHAVIOR OF YANG-MILLS 
GREEN FUNCTIONS 

Neglecting the trivial color factor 6 a b, the gluon and 
ghost propagators are parametrized as 



D^(p 2 ) = ( 5„ 



P^Pu\ Z(p 2 ) 



and 



D G ( P 2 ) 



G{P 2 ) 



P 



(1) 



(2) 



For very small momenta it is expected that the dressing 
functions behave power-like: 

Z IR {p 2 ) = c , 2 • (p 2 ) 2 -, G IR {p 2 ) = c 2 ,o • {p 2 )- K . (3) 

The exponents of the squared momentum are called in- 
frared exponents and depend only on one parameter, k. 
An approximate numerical value n « 0.5953 can be de- 
termined from the IR consistency relation between the 
ghost and gluon DSEs (6j [7J. Recently a similar value 
k w 0.57 was obtained in a strong coupling lattice anal- 
ysis [20]. In agreement with the Kugo-Ojima and the 
Gribov-Zwanziger scenarios the ghost propagator is di- 
vergent in the infrared, whereas the gluon propagator is 
slightly vanishing like (p 2 ) ' 2 . (NB: For an analysis of 
the K-dependence of the propagators see refs. j6] [19] .1 
The three-gluon vertex dependence was investigated in 
refs. |18l fl9] . In the uniform limit, where all momenta 
approach zero, also the dressing functions of vertices are 
given by general power laws. Each vertex with 2n ghosts 
and m gluons then has at least one dominant tensor that 
behaves as C2 njTO (n, {p 2 /p 2 }) • (jp 2 ) S2n - m with some pref- 
actor c<i n m , provided no non-trivial cancelations occur. 
&in.m is hereby the uniform IR exponent that determines 
the IR behavior in terms of one common momentum scale 
p\ 

The general expression for the (uniform) IR exponent 
in d dimensions is |18l [T9l [21] 



^2n,m = (n - m)K + (1 - n) 



-2 



(4) 



One possible way to obtain Eq. Q is via the non- 
renormalization of the ghost-gluon vertex in Landau 
gauge [22 , which is connected to a finite IR dressing 
of the vertex: 5 2 ,i — 0. This has been confirmed with 
DSE methods for d = 3,4 [23] and calculations on the 
lattice for d = 2, 3, 4 [M] [25] [26] [27] [28]. For higher ver- 
tex functions, e.g. the 1PI five-gluon vertex, a skeleton 
expansion was used in refs. |18l H"9l [2~T] . A feature of this 
expansion is that all orders possess the same IR expo- 
nent [18], which is true also for two and three dimensions 
|19l 121] . As a matter of fact the IR exponent is equal for 
all diagrams with the same type of bare vertex [19] . For 
the three-gluon vertex DSE diagrams containing a bare 
ghost-gluon vertex turn out to have an IR exponent of 



—3k. Those with a bare three- or four- gluon vertex have 
and ft, respectively. Another feature of this solution is the 
IR fixed-point of the running coupling |5l [T8l [T9l [29l [30] . 

In all these considerations it was assumed that ver- 
tex functions feature only singular behavior when all mo- 
menta vanish uniformly. Taking into account that ver- 
tex functions depend on several independent momenta, a 
more refined picture emerges that does not alter but ex- 
tend the previously described one [2]. The dependence 
on different momenta is encoded in the explicit depen- 
dence of the vertex dressing functions on all indepen- 
dent momenta. A general vertex function can then show 
power law scaling with different IR exponents 6ij de- 
pending on the momentum configuration pj (g 2 , • • • , (j, 2 ) 
that vanishes: 

T^-^( qi ,... ,q n ) = 



Such a non-trivial momentum dependence is also found in 
the analytic IR solution and is in perfect agreement with 
the power counting analysis performed in ref. [11], where 
all IR exponents for the primitively divergent Green func- 
tions were determined. For the three-point functions 
three new IR exponents for the case when one momen- 
tum becomes small compared to the other two were intro- 
duced. In the following we present explicit calculations 
that confirm the related previously obtained results of 
ref. [H]. 



Sgh 


5gl 








u 33 


°3g 


— K 


2/t 





— 3k 





min(| - 2k, 0) 


min(l - 2k, 0) 



TABLE I: The IR exponents for the ghost propagator, the 
gluon propagator, the ghost-gluon and the three-gluon ver- 
tices with all momenta going to zero, the ghost-gluon ver- 
tex with a ghost momentum going to zero, the ghost-gluon 
and the three-gluon vertices with a gluon momentum going 
to zero. If not denoted otherwise, we use k = 0.5953 ... in 
this article. 

The IR analysis of the three-gluon vertex DSE [18] 
identifies the diagrams containing ghosts as the IR lead- 
ing ones. The second ghost diagram of the three-gluon 
vertex DSE in Fig. [l]can be expanded using a skeleton 
expansion. Here the lowest diagram possesses two loops. 
The first order of the skeleton expansion is therefore the 
left ghost diagram in Fig. [l] called ghost triangle in the 
following. The same analysis for the ghost-gluon vertex 
DSE identifies the fourth diagram on the right hand side 
of Fig. [l] called ghost-ghost-gluon triangle, as the leading 
first order diagram. Both these diagrams are calculated 
in the subsequent sections. The complete DSEs as given 
in Fig. [l] can be derived by the algorithm explained in 
|31j . which is implemented in the Mathematica package 
DoDSE. 
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FIG. 1: The DSEs for the three-gluon and ghost-gluon vertices. They were derived using the algorithm explained in ref. |31| . 



Two notes on the validity of the results of the pre- 
sented calculations are in order: Because we do not solve 
a full system of DSEs but rather calculate single diagrams 
for low external momenta we only can provide qualita- 
tive solutions. Furthermore, we treat only the Yang-Mills 
sector. This is sufficient for an analysis of QCD in the 
quenched approximation [17 as well as in the unquenched 
case j!5j . 



III. THE GENERALIZED TRIANGLE 
INTEGRALS 

We will calculate the two diagrams given in Fig. [2] In 
the following it is understood that the DSEs have been 
renormalized by a suitable subtraction scheme [5] and 
no terms divergent in the ultraviolet are present. Thus 
all quantities are to be considered ultraviolet finite even 
without writing down explicitly the renormalization con- 



stants. Our presentation will mainly focus on the ghost 
triangle as the procedure is quite similar for the ghost- 
gluon vertex. Using bare ghost-gluon vertices, ig f abc 
with g M being the outgoing ghost momentum, but dressed 
propagators, the ghost triangle correction in the renor- 
malized three-gluon DSE yields 



d d q 



abc Nc c 2,o 



3 „3 



G((q + Pl ) 2 )G((q-p 2 ) 2 )G(q 2 ) 
(q + Pi) 2 (q~P2) 2 q 2 



(5) 



For small external momenta pi and p 2 the integral will 
be dominated by the singularities of the IR ghost propa- 
gators (Eq. |3|). It can be decomposed in one vector and 
two tensor integrals of rank two and one of rank three, 
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FIG. 2: Momentum routing for the ghost and the ghost-ghost-gluon triangles. 



= » /° 6C C ' { - P1^P2 P I»(P1,P2,P3; 1 + K, 1 + «, 1 + KJ d) + 

+ Pi li Iu P {pi,P2,P3> 1 + n,l + k,1 + n;d) + I^ P (pi,P2,P3; 1 + k, 1 + k, 1 + k; d)- 

- P2 p I^(pi,P2,P3; 1 + k,1 + k,1 + k; d))}, (6) 
where the tensor integrals are defined as 

I^...{puP2,P3^ lM] d) :=y (2^ ((g±p2)2) ^ ((g±pi) 2^ (9 2^- ( 7 ) 

The signs in the denominators depend on the momentum routing, but the solution for the scalar integrals is inde- 
pendent of them. The positive factor C2.0 in the ghost dressing function and the coupling g are to be considered as 
unknown constants because we have no scale to determine physical values for them. Therefore the integrals can only 
give information about the qualitative behavior of the three-gluon vertex as mentioned above. The factor —N c /2 
stems from the color algebra. In all calculations only the term in braces was considered. For the ghost-ghost-gluon 
triangle the same procedure was performed. 

Using a method due to Davydychev [32] we can decompose the tensor integrals into scalar ones, 

Pa; i*. «*; d) := / ^ {{q±n?)H{ l ± pi?)H<fY s ■ W 



The advantage of this method, described in more detail 
in App. [B] is that it involves only a minimal number 
of scalar integrals in higher dimensions. In contrast, a 
direct decomposition into four-dimensional integrals in- 
volves a vast number of different scalar integrals that 
mostly cancel mutually and yields equivalent results as 
we have checked in certain cases. 



The remaining scalar integrals involve denominators 
with non-integer exponents. The solution for such inte- 
grals was obtained both via a Mellin-Barnes representa- 
tion for occurring Gaussian hypergeometric functions [33] 
and within the Negative Dimensions Integration Method 
(NDIM) |34j . The result in general dimension, involving 
Appell's hypergeometric series F4, is given by 



l(pi,P2,p 3 ;vi,v2,v3;d) = —jipD^-^-^x 



5 



p2\ #-«2-"i / p 2x a r( | _ ^)r(-| + ^ + i^)r(-| + v 2 + u 3 ) 



y { \pi) \pi) ' n^n^n^) 

r~ i , d ^ d d p\ p\ 

x F 4 [d - v 2 - v 3 ~ v Yl - - v 2 ] 1 + - - v 2 - v\, 1 + - - v 2 - v 3 ; , -=■ 
2 2 2 

r(| - ^ - n)r(f - ^ - ^)r(-f + ^ 2 + ^ + ^ 
r(^ 3 )r(i/i)r(d - v 2 - v 3 - ^1) 

d d d p\ 

- + v 2 + v 3 + vi;l-- + v 2 + v u l-- + v 2 + v 3 ;^ i ,^ i 
2 2 2 pj. pi 



ply-" 2 - 1 ' 1 T(l - -v 2 - u 3 )T(-l + v 2 + pj) 

d d d p\ p\ 

- + ^1,-2^3; 1 + - +2/2 + ^1,1 - - - V2 - V 3 \ ~2, -2 

£ I I Pg Pg 

i-^-^ r( | _ ^ )r( d _ ^ 2 _ ^)r(_| + V2 + ^3) 



Piy r(z/ 2 )r(z/ 3 )r(d - f 2 - ^3 - ^1) 

x F 4 ( i^i, --2/ 3 ;i -2+^2 + ^1,1+ 2 - ^ - ^ 3 ;p,p ) } (9) 

\ ^3 ^3 / 



with Appells' series F4 defined as 



F 4 (a,6; c, d;x,y) 



x -» (a, m + to + n) x m y n 
~Ln (c,m)(d,n) ml n\ 



(10) 

The so-called Pochhammer symbol (a, n) is given in Eq. 
(All. When interested in the Euclidean momentum re- 



gion this result is not directly applicable, since the defin- 
ing series does not converge there. Its region of conver- 
gence, given by 



y\ <i, 



(11) 



is shown in Figj3]as the dark gray region. Via known ex- 
pressions [35, 36, 37 this series can be continued analyti- 
cally to the light gray areas which correspond to solutions 
for the two other possible ways of forming momentum ra- 
tios in eq. |9|. 

In contrast, using momentum and energy conservation 
for the three external momenta p\, p 2 and p 3 , the Eu- 
clidean region is defined by (a is the angle between p 2 
and p 3 ) 

P\ = (-P2 - P3 ) 2 =p 2 2 +pl + 2\ P2 \\p 3 \cos(a) (12) 

leading for x = p\jp\ and y = p\/p\ to the expression 

x + y - 2y/xy < 1 < x + y + 2^/xy. (13) 

This region is precisely the white part in Fig. [3] A 
series representation of Appell's function F 4 convergent 
in a part of the Euclidean region is given in Appendix [A] 
Employing symmetries and other, simpler continuation 
formulae one can cover the whole Euclidean region. 



0.25 0.5 0.75 1 1.25 1.5 1.75 2 



FIG. 3: The regions of convergence for different series repre- 
sentations of Appell's function F 4 . 



Recently it was claimed (38] [39] that Eq. (|9| is not 
correct and the number of Appell's series should be re- 
duced to three. However, the procedure performed to 
this end only amounts to a standard analytic continua- 
tion of two Appell's series to get only one. This is not 
per se forbidden, but the new Appell's series has a region 
of convergence (the upper light gray region in Fig. [3} 
that does not intersect that of the two remaining ones 
(the dark gray region), so that the total result does not 
converge. We therefore do not see, how such a reduction 
is of any help. 

An important feature of the three-point integral solu- 
tion are divergences at some kinematic points. Appell's 
function F4 may possess poles at (00, 00), (0, 1) and (1, 0), 
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but it can also be finite or vanishing depending on the 
parameters. Another reason for kinematic divergences 
are the prefactors in Eq. |9|, which diverge for certain 
values of the parameters in case of collinear momenta, 
i.e. one momentum is zero. Approaching these poles is 
equivalent to letting only one momentum vanish and is 
described by the non-uniform IR exponents when only a 
single momentum vanishes. We call these points asym- 
metric points. Since the final values of the IR exponents 
are a result of the intricate interplay between the diver- 
gences in the prefactors and of the Appell's functions, 
we extracted all values numerically and checked some of 
them analytically. As an example the quantitative result 
for the scalar integral I(p\,p 2 ,,p\ = 1; l + /t, 1 + K, 1 + k;4) 
is shown in Fig. [4] The analytic result is compared to 
a result from a numerical integration and shows good 
agreement. The general features of this function are that 
it rises monotonically from the symmetric point (1, 1) 
towards the boundary of the Euclidean region. At the 
boundary it has singularities at (1,0), (0, 1) and (oo,oo) 
and saddlepoints in between. One can clearly see these 
divergences when one momentum becomes small. 

A remark on the asymmetric point: Under certain cir- 
cumstances one can evaluate the integral exactly at the 
point pi = by use of the two-point integral, 



d d c 



\ <-> 



(q 2 )^ ((q±p) 2 )^ 



where the conditions 



Mi)r(| 



•A*a)r(/ii+A*2-|) 



r(/ii)r(/x 2 )r(d - hi- /i 2 ) 



(14) 



d/2 - m - ^2 < 0, c2/2-jUi>0, 



d/2 - a*2 > 

(15) 



have to be fulfilled for IR and UV convergence. One 
can directly show this expression to be equivalent to the 
three-point solution at the asymmetric point for \ — v\ — 
1/3 > 0. Most easily this can be seen for the case p\ — ► 
by setting it to zero in Eq. ([9]). Thereby the Appell's 
series become one-dimensional Gaussian hypergeometric 
series 2-P1 and two terms vanish if the exponents of p\ 
are greater than zero (d/2 — v\ — z/3 > 0). For p\ = 
p\ the Gaussian series have unit argument and converge 
only for 1 — Vy — 1/3 > 0. This means that for d > 2 
the latter condition is stronger. If it is fulfilled the two- 
point and three-point integrals coincide, which can be 
seen by setting \iy = v\ + V3 and /j, 2 = v 2 . The second 



condition in Eq. ( 15 I corresponds then to d/2 — v\ — 1/3 > 



0. The case when this condition is not true yields the 
kinematic singularities treated in this article and the two- 
point integral cannot be used for calculating this limit. 

As explained in Sec. [TlJ the power counting analysis 
only gives the IR exponent for the most divergent dress- 
ing function that is possible. Neither do specific results 
have to match the exponents nor can less divergent, finite 
or even vanishing dressing functions be excluded. Indeed 




FIG. 4: The scalar integral Eq. (|9| in the Euclidean momen- 
tum region at fixed momentum p 3 . Overlayed are both the 
analytic solution (black contour lines) and the result of a nu- 
merical integration (gray contour lines) for comparison. The 
slight deviations near the singularities result from the slow 
convergence of the numerical integration routine. 



the results presented below show that the combination of 
Appell's functions in the scalar integrals can lead to the 
realization of these additional possibilities. 

We would like to mention that in principle the inves- 
tigation of the four-gluon vertex is possible in an analo- 
gous way. But although NDIM enables one to derive a 
solution for the ghost-box integral, the situation is more 
involved: Without approximations or constraints on the 
external momenta fivefold hypergeometric sums emerge, 
which are technically much more difficult to handle than 
a double sum. Furthermore we expect that an analytic 
continuation will be necessary as well, but we know of 
none at the moment. Apart from these problems with the 
integrals the tensor structure is more complicated since 
the four-gluon vertex is constructed of a tensor of rank 
four in Lorentz space. Taking also into account that a so- 
lution in one region of convergence consists of more than 
four series 1 it is evident that in every step a multiple 
of the amount of work is necessary than for the three- 
point integral. Nevertheless calculations have been done 
in refs. [401 |4"T] for a given momentum configuration. 



1 The complete (naive) solution for the three-point integral as de- 
rived with NDIM comprises twelve series (three regions with four 
series each), whereas for the four-point integral this number is 
162. 
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IV. INFRARED SOLUTION OF THE 
THREE-GLUON VERTEX 

A perturbative solution for the three gluon vertex in 
arbitrary linear covariant gauge and dimension valid in 
the UV regime is given in [42 . Here we give a semi- 
perturbative result that involves fully dressed propaga- 
tors valid in the IR regime. With a general solution to 
the massless three-point integral (Eq. ^ together with 
the results from the analytic continuation presented in 
Appendix |A| and Davydychev's method (32] for calcu- 
lating tensor integrals (see Appendix |Bj at hand we are 
able to present here the result of the IR dominant part of 
the three-gluon and ghost-gluon vertices in the complete 
Euclidean momentum regime in the minimal tensor ba- 
sis provided by the employed computation method [32]. 
Furthermore, we can give the overlap of the vertex with 
arbitrary tensors and do so for the transverse part of the 
vertex for a comparison with lattice data. 

The results are obtained completely analytically and 
allow in principle to get the IR exponents directly. How- 
ever, due to the involved form of the integrals it is much 
faster to extract them by numerical fits. We did this for 
< K < 1 and d — 2, 3, 4 dimensions and present the re- 
sults in table |III| which also shows that the dependence 
on k is more involved as originally anticipated, since be- 
sides the additional term arising in two and three dimen- 
sions it is also possible that the dressing functions become 
constant. In the three-dimensional plots and Fig. [8] we 
used k = 0.5953 . . ., which is the best currently available 
value EE]. 

A. General Kinematic Dependence 

As explained in Appendix [B] since we only take into 
account the tree-level tensor of the ghost-gluon vertex, 
for the IR part of the three-gluon vertex only ten of the 
possible 14 tensors arise: 

10 

i=l 

(16) 

with dressing functions Ei and tensors r* as given in 
table [IT] The solution involving the scalar integrals as 
given in eq. (|oJ consists only of hypergeometric functions 
which can be computed, if converging series representa- 
tions can be found. Appendix [A] provides such represen- 
tations for the three-gluon vertex. The calculation of the 
series was performed with Mathematica 6 [43,. Conver- 
gence of the series is very good except at the boundaries 
of the Euclidean region. The free variables of the solu- 
tion are pf, p\ and p|, which are combined to x = p\jp\ 
and y — p\jp\- Due to the symmetry the full solution 
is contained in the unit square and is then parametrized 
by one unbounded scale and two bounded momentum 
ratios. To each point in the x-y-plane corresponds an 




2.0 



FIG. 5: For p§ = 1 the angle between pi and p2, which is given 
in general by cos a = (1 — a; — y)/(2y/x y). The boundaries 
represent the collinear case pi ~ P2 and the dot lying on the 
27r/3 contour line represents the symmetric configuration p\ = 

2 2 
P2=P3- 

angle between p\ and pi which is depicted in Fig. [5] to 
show the connection to another parametrization of the 
independent quantities employing two absolute values of 
momenta and the enclosed angle. A change of p\ corre- 
sponds to moving the kinematic point in the (x, y)-plane 
along radial lines, whereas scaling all momenta by the 
same factor does not change it. Thus uniform scaling 
Pi ~ p\ ~ p\ = P 2 with all ratios fixed and finite does 
not depend on the hypergeometric functions, where only 
x and y appear, but only on the prefactor in Eq. 
Taking into account the momenta from the bare ghost- 
gluon vertices and the anomalous scaling of the ghost 
propagators this gives in four dimensions 

(p 2 )!-^^ 2 - 1 ' 3 -> ( p 2)2-3(l+ K )+3i = (p2)i-3* ( 17 ) 

For the three-dimensional plots we fixed p\ = 1. This 
sets x and y equal to p\ and p\. In the following we will 
treat the four-dimensional case with k — 0.5953, except 
where denoted otherwise. 

Although the full solution of the three-gluon vertex is 
fully symmetric the choice of a tensor basis can hide this 
symmetry. The reason for choosing the solution of the in- 
tegrals with p\ in the denominator of the variables of the 
F4 function is that our basis depends only on p\ and p2- 
Having the same independent momenta for the basis and 
the variables of the Appell's function explicitly reveals 
some symmetries under exchange of the two momenta in 
the plots, e. g. the first and the second scalars corre- 
sponding to the tensors Pi^PuPiJpj and P2 li P2 v P2 p lp\- 
The sign is different because of the color factor. 

The results for the ten dressing functions are given in 
figs. [6] and [7] Although they only mildly change near 




the symmetric point, p\ = p\ — Vz-, there is a strong 
dependence in the region near the boundary correspond- 
ing to the collinear kinematics p\ ~ p2 and in particu- 
lar when approaching the singularities at the asymmetric 
points pi = 0. In contrast, previously used parametriza- 
tions of the vertex took only into account its scale de- 
pendence and therefore a constant function in the corre- 



sponding figures. To study the behavior near the asym- 
metric points in more detail we calculated the line start- 
ing at the symmetric point as shown in Fig. [8] It has 
the advantage of fewer terms necessary to be calculated, 
because y = 1, and we can get as far as p\jp\ — 1CP 8 . In- 
terestingly the power law behavior sets in for rather low 
ratios of p\/p\, but in the region below 10~ 3 the power 
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laws behavior can clearly be recognized. At each asym- 
metric point five dressing functions are diverging whereas 
the remaining ones simply become constant or even van- 
ish. 

We can get values for the exponents by calculating the 
approach to the asymmetric points. Fitting the expo- 



nents according to a power law and expressing the result 
in terms of k = 0.5953 gives the values reported in table 
[iljwith rather mild kinematic singularities with exponent 
1 — 2k. Note that the only contributions that really di- 
verge in the case of small p\ stem from S7T7 and EaoTio, 
since the other divergent scalars are rendered effectively 
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vanishing when multiplied by the corresponding tensors. 
However, as in Landau gauge the gluon propagator is 
purely transverse, the longitudinal parts of vertices do 
not contribute in DSEs [9]. The tensors of Ej as well as 
of the part of E\q that is not suppressed by pi are longi- 
tudinal with respect to p\ and p$ respectively. Therefore 
we do not find divergent transverse parts of this diagram 
that contribute in any DSE. This is in agreement with our 
previous power counting analysis [14J: If the external legs 
of the ghost triangle in Fig. [2] are projected transversely 
the canonical dimensions of all three ghost-gluon vertices 
scale with the soft momentum, so that the expected IR 
exponent is 

<t = \ + \ + 2(V - 1) + 2(S gg + \) = | - 2k (18) 

in contrast to the non-projected case, where we had 
1 — 2k. This means that in the transversely projected 
case the bare vertex is leading in the soft gluon limit, 
^ig t — 0- However, if the vertex is contracted with some 
tensor, cancelations may occur so that the ghost-triangle 
is dominant. We show in subsec. IIV Bl that this is in- 
deed the case and it should be possible to observe the 
kinematic divergence of Eiq on the lattice due to an addi- 
tional suppression of the tree-level part at the asymmetric 
points, when contracted with a certain tensor. Since we 
do not calculate all 14 tensors, it may be possible that 
adding them might change this special role of Eiq. We 
also neglected the one-loop ghost term (third diagram on 
the right-hand side of Fig. [TJ in the three-gluon DSE 
that is a leading term in the conformal limit, but one 
order higher in the skeleton expansion. 

In table [II] we give the IR exponents of all ten tensor 
components for the physically relevant case d = 4 and 
k = 0.5953 |6j 13- Interestingly the momentum depen- 
dence is not independent of k, as we checked with explicit 
calculations for values in the range [0,1]. The result can 
most easily be presented in a closed form when using 
unnormalized dressing functions: Then one can see that 
when the IR exponent is zero it does not become larger 
for lower values of k but stays zero so that the dressing 
functions are constant. Due to different normalizations 
of the tensors this happens at different values. For ex- 
ample in four dimensions Eio becomes finite for k < 1/2, 
whereas Ex is still divergent. For this reason we present 
the complete result valid for d = 2,3,4 and < n < 1 
in terms of unnormalized tensors in table IIIII The addi- 
tional term (d — 4)/2 accounts for the dependence on the 
dimension. The extension to two and three dimensions 
might prove useful for instance for comparison with cal- 
culations on the lattice, where often lower-dimensional 
studies are performed, e. g. |27| [28] . 

Let us now discuss the qualitative features of these 
results that should be independent of the considered 
truncation and their impact on physical quantities. 
The three-gluon vertex enters for instance in the gluon 
loop contribution in the DSE for the gluon propagator. 
For soft external momenta in the deep IR regime but 



hard loop momenta the kinematic divergence is directly 
probed. Despite their much larger support in the loop 
integral (two momenta can be arbitrary and only the 
third one has to become small) the mild soft-gluon di- 
vergence is too small that it could be relevant and is IR 
suppressed compared to the leading ghost loop as has 
been discussed in [T3]. In contrast for finite momenta 
of the order of hadronic scales the gluon loop also be- 
comes important. Here the kinematic divergences are 
not directly probed. However, the results show that they 
also influence the behavior for finite momenta and in- 
duce a strong dependence on the kinematics in many 
dressing functions. The generic feature in this case is 
that even for finite momenta their absolute values in the 
vicinity of the boundary of the Euclidean region, corre- 
sponding to the collinear kinematics, and in particular 
towards the asymmetric points are considerably larger 
than around the symmetric point. Of course, at hadronic 
scales our results are not complete due to neglected glu- 
onic contributions to the three-gluon vertex that may be 
even of comparable size in this regime. Yet, the found 
pronounced kinematic dependence is not expected to be 
completely washed out in a more refined analysis. Corre- 
spondingly this kinematic dependence should be relevant 
for the quantitative behavior of the gluon propagator in 
the intermediate regime, which in turn has a dominant 
impact for the properties of bound hadronic states within 
Bethe-Salpeter and Fadeev approaches, see e. g. [8]. Sim- 
ilar remarks hold for the quantitative influence on the 
quark-gluon vertex and the corresponding impact on the 
recently proposed confinement mechanism |15l I17J . 



r ~ -Sl 




FIG. 8: Approaching the asymmetric point (0, 1), i. e. p\ — > 0. 
Five scalars diverge, three are constant and two vanish. From 
top to bottom at the left side the scalars are: Eg, Eio, E$, Et, 
Ei, E2, E4, Eq, Eg, E%. The corresponding case for vanishing 
p\ can be inferred from the Bose symmetry of the three-gluon 
vertex. 



B. Three-Gluon Vertex on the Lattice 

Since this article is the first that gives qualitative 
numerical results for kinematic IR singularities it is of 
course of interest to investigate the possibilities of other 
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Tensor 


pi 


soft 


P2 


soft 


r i f ^P = 




1 - 


- 2k 





T2^ p = 


P2 p .P2„P2 p /pl 





1 - 


- 2k 


T3 fl „ p = 


(Plv.PlvP2 p + Pl fI P2„Pl p + P2 p Pl„Pl p )/p? 


2 - 


- 2k 





T W = 


(Pl„P2„P2 p +P2„Pl„P2 p +P2 p P2„Pl p )/pl 





2 - 


- 2k 


T5 fl „ p = 


g^Pi p + g^Piv + g^pPi,, 


1 - 


- 2k 





T&^„ = 


g^P2„ + 3mpP2„ + SWP2 P 





1 - 


- 2k 


T W = 


Pl M Pl„P2 p /Pl 


1 - 


- 2k 


1 - 


- 2k 


= 


Pl p P2„P2 p /pl 


1 - 


- 2k 


1 - 


- 2k 


T"W = 


(Pl M P2„(P2 ~Pl)p + (P2 -pi) /J Pl„P2 p )/(piP2) 


3/2 


- 2k 


3/2 


- 2k 


no M „ p = 


= gvpPi» - g^P2 p 


1 - 


- 2k 


1 - 


- 2k 



TABLE II: The IR exponents in terms of k for the normalized scalar functions corresponding to the indicated ten sors of the 
ghost triangle when p\ or p\ become soft in four dimensions. Since the dependence on k is not trivial (cf. table III I, these 
results are only valid for k > 1/2. 



Tensor 


pi soft 


P2 Soft 


= Pi p Pi„Pi p 


-2k+^ 





T~2 pvp = P2 P P2„P2 P 





-2k +^ 


=Pl p Pl„P2 p +Pl t ,P2„Pl p + P2 p Pl v Pl p 


min(0, 1-2k+^) 





74 p „ p =Pl p P2„P2 p +P2 p Pl„P2 p +P2^P2„Pl p 





min(0, 1-2k + ^) 


^5 P „ P = SWPi p + SWPi„ +3^pPl p 


min(0, 1-2k+^) 





fe^p = g^P2 p + g^ P P2 v + g vp p2 p 





min(0, l - 2k + <±±) 


?W = Pl p Pl„P2 p 


-2k +^ 


min(0, 1-2k + <*= i ) 


= Pl p P2„P2 p 


min(0,l - 2k + ^) 


-2k +^ 


TB*lup =PI P P2„{P2 ~Pl)p+ (P2 -Pl) M Pl„P2 p 


min(0, 1-2k + 


min(0, 1-2k+ *=i) 


T10^ P = g vp p\ p - g^P2 p 


min(0,l- 2k + ^) 


min(0, 1-2k+ ^) 



TABLE III: The IR exponents for the unnormalized scalar functions in terms of k corresponding to the indicated tensors of 
the ghost triangle when p\ or p\ become soft. In two and three dimensions additional contributions arise. The choice of 
unnormalized tensors is only for an easier and more intuitive representation of the results. For comparison we give the values 
of k in two, three and four dimensions that are currently believed to be the most reliable ones: 0.2, 0.39 . . . and 0.59 . . . 
respectively |6l 171144]. 



methods to confirm their existence. One natural possibil- 
ity are lattice Monte-Carlo simulations. In this subsec- 
tion we want to compare our results with recent lattice 
studies and investigate, if there exists a possibility to ob- 
serve kinematic singularities on the lattice. 



In general a comparison can only be done for the four 
transverse components of the three-gluon vertex, i.e. the 
tensors t\, T2, tj and r§ will give no contributions to 
quantities calculated on the lattice. Furthermore, a di- 
vergence would not manifest itself in the same way as in 
the continuum, because on a finite lattice results are nat- 
urally always finite. Thus, calculations can be done di- 
rectly at an asymmetric point and singularities can only 
be seen in the volume dependence of the investigated 
quantity. Although there have been calculations directly 
at an asymmetric point |25l l27l 128] . they have not found 
any indication of a singularity. The quantity used in refs. 
|251 [28] was a projection of the full three-gluon vertex to 



the tree-level expression normalized appropriately: 

q = r|Sp {pi )D^ a (pi)D„p(P2) Dfn p3 )r a/ 3 7 {Pi ) ^ 

r|S P fa )D fia (p 1 )D L ,p(p 2 ) d p1 (p 3 )r^ 7 (pi ) 

Gluon propagators are denoted by D and the dressed and 
bare three-gluon vertices by V and r' ', respectively. Our 
results indicate that no divergent contribution from the 
ghost triangle can be expected for p\ — p\ and p\ —> 
and there is no contradiction between lattice studies and 
our results from the DSE approach as far as the three- 
gluon vertex is concerned. It can be shown analytically 
that the only surviving scalar in this limit is Eq, which 
is finite at p\ = 0. The corresponding numerical result is 
shown in Fig. [9] 

Projections with other tensors may yield different re- 
sults. To investigate this possibility we resort to the Ball- 
Chiu basis [45], which explicitly provides four transverse 
tensors. These are a natural basis for future studies of 
the full three-gluon vertex on the lattice, where only the 
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transverse part of the vertex is accessible: 

F jlu P (PUP2,P3) = (<Wl -P2 -Pl„P2jB% 
F lvp(P2,P3,Pl) = (8vpP2 -P3 ~P2 p P3 v )Bl, 
F^p(P3,Pl,P2) = (Sf,pP3 - Pi -P3„Pi p )Bl, 
H^p(p!,P2,P3) = (SpuBp - KpB\ - 5p,pBl- 
+ Pl p Pl„PZ„ ~PUP2 P P3 P ), 



(20) 
(21) 
(22) 

(23) 



with £?* = tijkP 3p Pk ■ Pi- The tensor H^ vp is the only 
one that depends on pi in every term, so that it scales 
like (pf) 1 / 2 in the corresponding limit. The calculated 
contractions are: 



1 



F l {pi,pi,Pt) = -^F^ p (p 1 ,p 2 ,p 3 )P fla P l/0 P pi T 



Pi 
1 



a/37 5 



H {P 1 ,P 2 ,P3) = -2 H ^p{Pl^P2,P3)PpaP V fiPp^ a 



Pi 



01 > 



(24) 



where P M „ is the standard transverse projector and 
the ghost triangle part of the three-gluon vertex. The 
results are given in Fig. 10 The normalization 1/p 2 was 
chosen to show the qualitative different behavior of the 
F l and the H quantities. Analytic considerations show 
that directly at the asymmetric point, i. e. p\ = and 
P2 = the corresponding values are 

F x {Q,plpl) = 3(pl) 2 E 6 (Q,plpl), (25) 
F 2 (0,p 2 ,p 2 ) = F 3 (0,plp 2 2 ) = E 6 (p 2 2 ) 2 (0,pIp 2 2 ), (26) 
H(0,p 2 ,p 2 ) = 2E w p 2 (0,p 2 ,p 2 ). (27) 

Thus the F % are finite and H diverges like (p 2 ) 1_2fc . Nor- 
mally one would expect that the tree-level gives the dom- 
inant contribution, since it does not scale in the soft 
gluon limit, whereas the transverse part of Eiq is sup- 
pressed. However, it turns out that in the case p\ = p\ 
the contraction of the tree-level with the £f-tensor yields 
an additional suppression, i. e. although one expects 
that HPPPTW scales like (p 2 ) 1/2 , because of the p Y - 
dependence of the ff-tensor, it goes like (p 2 ) 1 - With 
normalization as given above this yields an IR finite con- 
tribution from the tree-level 1^°) and the ghost-triangle 
r A is the leading part in the contraction of the full ver- 
tex r = r' ) + r A with the 77-tensor giving a power law 
with a non-integer exponent. Thus it should be possible 
to verify the existence of kinematical singularities also on 
the lattice by studying the full vertex. However, we want 
to stress that only for p\ — p\ the leading part obeys a 
power law with an exponent that depends on n due to 
an additional suppression of the ff-tensor, whereas oth- 
erwise the tree- level is dominant. 



V. TENSOR COMPONENTS OF THE 
GHOST-GLUON VERTEX IN THE INFRARED 

The ghost-gluon vertex can in principle be treated 
along the same lines as the three-gluon vertex. One might 




FIG. 9: The normalized contraction of the ghost triangle of 
the three-gluon vertex with three propagators and a tree-level 
vertex as given in Eq. (19 1. The value clearly approaches a 
constant for p\ — > 0. 




FIG. 10: The contractions of the ghost triangle with t he t rans- 
verse tensors of the Ball-Chiu basis as given in Eq. (24 1 and 



the contraction of the tree-level tensor with the H tensor along 
the line from the symmetric to the asymmetric point p\ = 0. 
The values for F2 and F3 coincide. 



even suspect that the actual calculation is simpler, be- 
cause the full vertex only has two possible tensors, 

^^Pgh\Pgh\Pgl) = 

= iAiip^f^pDp^ + iB((p^) 2 ,p 2 gh ,p 2 gl )p gltt) 

(28) 

where Pgh/Pg 1 ^ is the momentum of the outgo- 
ing/incoming ghost and p g i that of the gluon. However, 
there are some additional complications that are due to 
the uniform IR exponent of the ghost-gluon vertex being 
0. Thereby the combinations of propagator IR exponents 
can combine to in the Gamma functions of Eq. 
This complicates the analytic continuation, because the 
limit has to be taken carefully in which case the diver- 
gent Gamma functions cancel each other to leave a finite 
result. The corresponding formulae that have to be used 
are given in appendix [A"] 

From the two existing DSEs for the ghost-gluon vertex 
we choose the seemingly more complicated one given in 
Fig. [T] where the bare vertices are attached to the ex- 
ternal gluon leg. In this case the leading diagrams are 
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Dressing function 


Pgl 


-> 


Pgh -» 


A A 


min(0, 2 - 




1 


B A 


min(0, 1 - 




1 


tree-level 








TABLE IV: The IR exponents of the ghost-ghost-gluon trian- 
gle and the tree-level when p g i or p g h become soft. 



Dress, fct. 


X tensor 


Pgl 


-> 


Pgh — ► 


Ap g h u 


min(0, 2 - 


-2 K +^) 


1/2 


B P9i^ 


min(|, | - 


-2/t+^) 


1 



TABLE V: The overall IR exponents of the ghost-gluon ver- 
tex, i.e. ghost-ghost-gluon triangle and the bare vertex in- 
cluding the canonical momenta. The dominant part of the 
transverse dressing function A is then given by the bare ver- 
tex. 



the ghost-ghost-gluon triangle and the ghost loop (dia- 
grams 4 and 5) as well as the bare vertex itself. The 
first order of the skeleton expansion consists of the first 
diagram only, since the ghost-ghost scattering kernel in 
the second is one-particle irreducible. The alternative 
version of the DSE features a leading diagram involving 
a dressed three-gluon vertex, whereas here only dressed 
ghost-gluon vertices appear and we employ again bare 
vertices. We denote the contributions to the dressing 
functions A and B from the ghost-ghost-gluon triangle 
by A A and B A respectively. 

The results for the ghost-ghost-gluon triangle in the 
Euclidean region are shown in Fig. [IT] The approach to 
the asymmetric points is plotted in Fig. [12] Here we have 
to investigate both limits of the gluon and the ghost leg 
momentum going to zero because there is no Bose sym- 
metry as in the case of the three-gluon vertex. Numeri- 
cal values for the exponents are given in table [TV] for the 
ghost-ghost-gluon triangle and in table [V] for the overall 
vertex. In four dimensions only the longitudinal dressing 
function B exposes a divergence when the gluon momen- 
tum goes to zero. Therefore the existence of kinematical 
divergences in the ghost-gluon vertex cannot be investi- 
gated on the lattice. Furthermore the divergent dressing 
function is multiplied with a tensor that also scales with 
the gluon momentum, rendering this contribution to the 
vertex IR vanishing in this limit. Our results agree with 
the statement of Taylor [22 that for vanishing ghost mo- 
mentum the ghost-gluon vertex becomes bare, i. e. the 
IR exponent is 0, because the loop corrections in the DSE 
are indeed subleading compared to the bare vertex. The 
dependence on k is similar to that of the ghost trian- 
gle, i. e. the IR exponent can also become for values 
of k below a certain value determined by the dimension. 
Whereas in four dimensions the transverse dressing func- 
tion A A is always finite for values of n between and 1 
it turns out that there is indeed a dependence on k that 
is revealed in lower dimensions. However, the relevant 
values of k are well below this regions, see table |V| 




FIG. 11: The dressing function contributions from the ghost- 
ghost-gluon triangle corresponding to the tensors p 9 ?i (ghost 
leg) and p gl ^ (gluon leg). 



VI. SUMMARY 



We presented an analytic solution for the IR leading 
parts of the three-gluon and ghost-gluon vertices, which 
are given by the ghost triangle and the ghost-ghost-gluon 
triangle. We found the expected uniform IR behavior, 
i.e. an IR exponent of —3k for the three-gluon vertex 
and for the ghost-gluon vertex, but also confirmed the 
existence of kinematic singularities in the case that only 
one momentum goes to zero as predicted in ref. [13]. As 
detailed above the non-trivial momentum dependence in 
the Euclidean region may have a quantitative impact on 
other Green functions. 
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For our calculations we combined known results for the 
massless three-point integral and analytic continuations 
of Appell's function F4. Using bare ghost-gluon vertices 
we found that for the three-gluon vertex five of ten dress- 
ing functions are divergent for one momentum going to 
zero. Although in this limit all transverse parts are addi- 
tionally suppressed, it should be possible to confirm the 
existence of kinematic divergencies on the lattice, since 
the otherwise leading tree-level tensor is additionally sup- 
pressed at the asymmetric point. 

Also for the ghost-gluon vertex there is a divergent 
dressing function that is suppressed by the corresponding 
tensor. Thus the vertex itself shows no kinematic singu- 
larities. We determined the dependence of the dressing 
functions on k in the region < k < 1 for two, three and 
four dimensions. 



1. Preliminaries 

The derivation is simplified by use of the Pochhammer 
symbol, defined as 



(a,n) 



T(a + n) 



(Al) 



Since it will be used extensively in the following deriva- 
tion, we give here some standard relations. Simply by 
inserting one can show 



(a, m + n) — (a, m)(a + m,n). 



(A2) 



From Legendre's duplication formula for the Gamma 
function a similar relation for the Pochhammer symbol 
follows: 
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APPENDIX A: ANALYTIC CONTINUATION OF 
APPELL'S FUNCTION F 4 INTO THE 
EUCLIDEAN MOMENTUM REGION 

In this appendix we derive the analytically contin- 
ued expression of Appell's function F4 that was used in 
the calculations for three-point integrals. The original 
derivation was performed in ref. [37] • We noted that this 
reference contained a few typos which we will correct 
in the following. In addition we discuss an appropriate 
reformulation for the treatment of special cases of the 
arising series which are necessary for a numerical imple- 
mentation. I 



(a, 26) = 2 zb (a/2, 6)(l/2 + a/2, b). (A3) 

An important formula is the analytic continuation of the 
Pochhammer symbol. It is only valid for integer values 
of n: 



(o,-n) 



(-1)" 
(1 — a, n) 



(A4) 



The Gaussian hypergeometric series and generalized hy- 
pergeometric series of one variable will appear, which are 
in general defined as 



pF q (ai,...,a p ;bi,...,b q ;z) = p F q L 1 ''"',*" z 

\h,---,b q ; 



E 



(ai,n) ...(a p ,n) z n 
n) . . . (b q ,n) n\ ' 



p,qeN . (A5) 



If p = 2 and q = 1 this is the Gaussian hypergeometric 
series of dimension one, where by dimension we mean the 
number of variables (here z) . The and b. t are referred 
to as parameters. 
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2. Analytic Continuation 



The first step is to rewrite the original definition of F4, Eq. (10 1, in such a way that it contains the Gaussian 
hypergeometric series of one variable: 



7-, / . , \ (a,m)(b,m) x n 

i<4 : {a,b;c,d;x,y) = — , \ ZT 2^i( a + m , ° + m\ a; y). 



m,n— 



(c, to) to! 



(A6) 



Now we use an analytic continuation of 2-F1 ((15.3.6) in |46]). 



„ , , \ T(c)T(c — a — b) „ . , 
2 F 1 (a, b; c; y) = — -4- 2 Fi(a, b; a + b - c + 1; 1 - y)+ 



r(c-o)r(c-6) 



r(c)r(a + fr- c ) 
r(a)r(6) 



(1 - 2/) c - Q - b 2 Ji(c - a, c - 6; c - a - b + 1; 1 - y). 



(A7) 



We insert eq. (A7| into eq. (A6| and treat the two appearing terms separately. The first one will lead to a series 
which we call G according to the conventions of ref. [37]. It is the easier one of the two emerging series: 

x m (a, m)(b, to) T(d)T(d — a — b — 2m) (1 — y) n (a + to, n)(b + m, n) 



E 



to! (c, m) T(d — a — m)T(d — b — to) 



E 



n\ (a + b + 2m-d+l) 



T(d)T(d~a-b) 



T(d-a)T(d-b) ^ to! n! 



E 



x m (1 — y) n (a, to + n)(b, to + n) (1 — d + a, to)(1 — d + b, to) 



(c, to) 



(1 - d + a + b, 2m + n) 



T{d)T{d-a-b) 
T{d- a)T(d-b) 



G(a, b, 1 — d + a, 1 — d + b; 1 — d + a + b, c; x, 1 — y). 



We have used Eqs. ( |A4| and ( Al I and the G series is defined as 

x m y n (a, to + n) (b, m + n) (c, to) (d, to) 



G(a,b,c,d;e,f;x,y) := ^ 



(e,2m + n)(/, to) 



(A8) 



(A9) 



The second part is more intricate. Employing Eqs. (A2|, (A4l and (Al I we get 

y x m (a, m)(6, to) r(d)I> + & - d + 2m) _ d _ a _ 6 _ 2ro>< 
•^-^m! (c, to) r(a + TO)r(6+m) 



E 



(1 — y)™ (d — a — to, n)(d — b — to, n) 
n! (1 + d — a — b — 2to, n) 



T(d)T(a + b-d) 

r(o)r(6) 



(1-2/) 



d—a—b „ 



S a 



^ V(i-y) 2 



(y - 1)" (a + b - d, 2m - 71) (1 + a - d, m)(l + 6 - d, to) 
n! (1 + a — d, m — n)(l + b — d, m — n)(c, m) 



(A10) 



In this form the series is still not convergent in the Euclidean region and we need to perform another analytic 
continuation. For this we rewrite eq. (AlOl into a form, which contains the generalized hypergeometric series 4F3: 

T(d)T(a + b-d) n „ d _ a _ b ~ (1 - „)» (d - a, n)(d - b, n) ^ 



<i-y)° 



E 



T(a)T(b) v ^ n\ (l-a-b+d,n) 



x 4 F 3 



b d n 1 I a i 6 

2 ' 2 + 2 + 2 2 

c, 1 + a — d — n, 1 + 6 — d — n; 



2 ' 2 2 2 ' 2 ^ 2 I 2 2 2 ' ^ I ^ ^7 1 I & C?j 



4:r; 



(i-y) 2 



(All) 



Eqs. (A2l, (A3 1 and (A4l were used. Writing the generalized hypergeometric series 4F3 as a Meijer-G function G^ n 
(eq. (rjon page 215 in |35j), 



iF 3 (a,b,c, d;e,f,g;x) 



r(e)r(/)r( g ) 

T(a)T(b)T(c)T(d) 



G 



4,1 



l,e,f,g\ 
a,b,c,d I ' 



(A12) 
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we can employ its analytic continuation (eq. (5) on p. 208 in |35j): 

i m 

n . = r(b j -b h )uu r ( i+b h- a ^ 



p,q 



a 1; . . . , a p 



E 



bi,...,b g j ^ nUn+i r d + h - bj) rru +1 r(aj - b h ) 

1 + bh - ai, . . . , 1 + bh - a. 



x bh x 



\1 + bh - h, ...,*,..., 1 + 6 ft - 6 g ; 



— m—n _ 



p < g Vp= g A |x| < 1. 



(A13) 



The primed product excludes the expression T(0) (when b = j), and the * stands for the one expression taken out as 



argument of p F q _i. Inserting eq. (A13I into eq. ( A12 I shows that two of the four expected 4.F3 vanish due to the 
appearance of T(—n) in the denominator which yields because n is an integer. In the two functions left half-integer 
values of n appear: 



T(d)T(a + b-d) 



T(a)T(b) 



(1-2/) 



d—a—b . 



r(c)r(|)(-4a ; )f- 



b _ d 
2 2 



TV E (-*) 



m\n 



Ax 



2' ra,n—0 

k _ Sl _ <L H _ mVS -I- - — - m — — 1 
2 2 2 ' 2 A 2^2 2 ' 2 / 



(r( c -|-|-|)r(i 

(l + f-|-|,f-m)(l 
(| ) m)(l + l-|-|,-f-m)(l-f + |-| ! -f-m)(c-f-| + |,f-m) 
(l-y)r(c)r(-|)(-4a;)^-f-S-^ ~ „ (-1)"+- ^(1-^ 

-p/„ _ a _ b , d _ l\T(a i b _ d\ Z^i ^ ' 
X V L 2 2^2 2' , l -2" r 2 2> m,n=0 



i.\n\ 



Ax 



(1 



b d n 

2 2 ' 2 



)(l+|-f-f ! f-m)(| + f + |-|, m -a) 



(f,m)(A-±a+f-f,-f -m)(l-f + §-f,-f -m)(c-f -§ + §-±,f -m) 



(A14) 



This expression could already be used in calculations. Nonetheless it proves convenient to rewrite the sums for 
easier implementation in a framework for treating hypergeometric series. Therefore the half-integer values have to be 
transformed into integer ones. This can be achieved by making the replacements ^ — » n and n — > n + | for the even 
and the odd parts respectively. Thereby one ends up with four different terms in which a new series called K can be 
identified: 



K(a, b,c, d;e,f,g,h;x,y) 



(a, m + n)(b 7 m + n){c,m — n)(d,m — n) x m y n 
^— ^ (e,m — n)(f,m — n)(g,m)(h,n) ml nl 



(A15) 



The result for the analytically continued series, called L, is 



L(a,b;c,d; — , -) := 

T(d)T(a + b-d) 1 d_ s _ § 

r(o)r(6) j 



i 



r(| 



f + i)r(c-f-| + f) 



/ 6 a _i_ d a 

x X 2 2 "I" 2' 2 



6 I d a | b d a r b 

2 i ~2'2 i ~2 2 1 2^2" 

6 a I d a b i rf 1_. 

2 2~ t_ 2'2 2~ t_ 2'2'2' 



+ !; (! — 2/) 2 x 
Ax ' 4 



(d+a-6-l)(d-a + 6-l)(-a;)3 
2(d - a - 6 + l)r(| + | - |)r(c- | - I + | + |) 



a I c/ 
2 f 2 



1 a _i_ b a! 

2 ' 2 2~ h 2~ h 2'2~ h 2 2 

a i d 1 tt ^ _i_ ^ 1 1 3 . 

2 2 + 2 2'2 2 + 2 2 ' 2 ' 2 ' 



1 a I b 

2'2 f 2 



4z ' 4 



b 1 



r( fl) rW ^ W-^d-y)^).-.-.-^ 
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X \r(§ + |-|)r( c -f-| + |-i) x 

(fa a | d I 1 a b | <i | 1 a i fa r 1 _i_ fa r d i 3 . f 1 7*^2 ™\ 

2 2 i ~2 i ~2'2 2'2'2'2'2 2 + 2 ' 2 + 2 6 2 + 2 > ^ y > _ | 1 

fa a I d I 1 a fa _r_ d 1 1 3 j^. 4™ ' 4 1 

2 2 ^ 2 ' 2 ' 2 2 ^ 2 ^ 2 ' 2 ' 2 ' / 

{d + a - b){d - a + 

+ 2(d - a - 6 + i)r(| + § - I - |)r( c - | _ I + |) x 

(fa _ a 1 d 1 1 a _ fa 1 <i 1 1 a , fa _ d ai6_„_d,-i.('-|_ 7( s i2 T \ | 

2 2 + 2 +1 '2 2 i ~2 +1 '2 i ~2 2> 2 + 2 C 2 + * > Z \ \ ( A 1 fil 

fa a, da fa , d 3 3 . 4x '41' ^ ' 

2 2 + 2> 2 2 + 2' 2> 2' ^ X V J 

From Eq. (All I to Eq. (Al6l Eqs. (A2l to (All were used extensively. Eq. (A16) is quite a lengthy expression and 
therefore Exton [37] abbreviated it with L(a,b,c,d; (1 — y) 2 /{4x), x/A) [37 . In contrast to the above defined G and 
K series, the combinations of x and y are fixed in the definition of L. This is caused by the appearance of x and y in 
the prefactors to the K series. L differs from the result given in |37J : 

• The second and fourth term have an additional factor 4. 

• The third argument of the second K series has an additional — A. 

• The second {—Ax) has an additional — A in the exponent. 
The final result for the Appell function is then 

F 4 {a, b; c, d; x, y) = ^4^ -4= " — 7V G ( a ^ fo > 1-d+a, l-d+b;l-d+a + b, c;x,l- y)+ 
1 (a — ajl (a — oj 

Ll^^cd;^ 1 ^ 2 ,^. (A17) 

I 



3. Regions of Convergence 

Now we determine the regions of convergence. For K 
we can employ the so-called method of cancellation of 
parameters. It can be shown |47] that the region of con- 
vergence does not depend on the parameters as long as 
these are not exceptional in the sense that the series is 
undefined, terminates or reduces to a sum of hyperge- 
ometric series of lower dimension. This means one can 
choose the parameters such that some Pochhammer sym- 
bols cancel each other. Choosing c — e and d — f the 
K series reduces to the standard F4 series for which the 
region of convergence is given in Eq. ( |Tl| . Inserting the 
occurring values for the variables the region of conver- 
gence for K, depicted in Fig. 13 (right), is 




< 1. 



(A18) 



For the G series Horn's theorem on the convergence of 
double hypergeometric series |48] has to be used. For a 
short overview see refs. |35l Wf\ . The resulting region of 
convergence is given in Fig. 13 (left). One can see that 
it contains the region of convergence of the K series and 
so the total region of convergence for Eq. (A17I is given 



by Eq. (A18I. There exist further analytic continuations 



region. They are given in the appendix of [37J and corre- 
spond to the analogous analytic continuations of the two 
additional series representations mentioned above (light 
gray regions in Fig. In principle they are not neces- 
sary here, because we can always choose the ratios of the 
momenta such that the area in the rectangle defined by 
(0,0) and (1, 1) is sufficient. Nonetheless we used several 
series in the numerical calculations because the conver- 
gence can be quite slow at the boundaries of the regions 
of convergence. One series used not given in [37 is the 
one obtained by exchanging c and d simultaneously with 
x and y. This is allowed because of the initial symmetry 
of F4 under this exchange. One should note that F 4 is 
not symmetric under the exchange of only c and d or x 
and y. The resulting region of convergence is Eq. (Al8l 
with x and y exchanged. 



4. Implementation of the Series 



General case: In the form of Eq. (Al6l the series 



representation is not suitable for a numerical evaluation, 
because it is possible that terms such as 



with which one can cover the whole Euclidean momentum 



(0,m— n) 

Depending on the values of m and n this can be or a 
finite value. In order not to have to treat all the possi- 
ble cases of peculiar combinations on their own a slightly 
different form for L was used. This is possible because 
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some parameters of the K series coincide. For the first 
and the third K series these are the first /second and the 
fifth/sixth parameters. Exploiting this K can be rewrit- 
ten as 



FIG. 13: The left and right pictures shows the regions of con- 
vergence for the G and K series respectively for the variables 
as appearing in Eq. ( A17 1, The region of convergence for the 



K series is the smaller one and thus K is the limiting series 
on the region of convergence. 



Ki,z{a, b, c, d; a, b, g, h; x, y) = 

r(a + 77i + n)T(b + m + n)(c,m — n)(d, m — n) x m y 

m,n— 



r(a + m — n)T(b + m — n)(g, m)(h, n) ml 77' 



(A19) 



In the second and fourth K series the fifth/sixth parameter differs from the first/second by one. As the parameters 
also appear in the prefactor we can combine all this to 



2(a — 1) 2(6 — 1) .^2,4(0, b,c,d;a — 1,6— 1, g, h; x, y) = 

r(a + 777 + n)T(b + m + n)(c, m — n)(d, m — n) x m y n 
r(a — 1 + 777 — 77)r(6 — 1 + 777 — n)(<7, m)(h, n) ml nl 



4 E 

m.n— 



(A20) 



Using Eqs. (Al9l and (A20l we did not encounter the 



problems appearing for the original K series. 

Special case: Normally Appell's function 
^4(0, 6; c, d; x,y) is unity if one of the two first pa- 
rameters vanishes. However, if it is multiplied by T(b), 
and b vanishes, the l/r(6) in Appell's function is can- 
celled and the formerly vanishing terms become finite, 
whereas the first summand is r(0). When calculating 



the ghost-gluon vertex exactly this happens for four 
integrals: The finite contributions become important 
and the T(0) terms cancel. To get the finite parts one 
has to repeat parts of the derivation from above. The 
bottom line is that some single summands are no longer 
part of the series and the analytic continuation has to be 
done very carefully as it only works for complete series. 
Right at the beginning we have instead of Eq. ( A6 1 



r(6)F 4 (a, b; c, d; x, y)\ b=a = T(0) + ]T ~ + E ^^[^ — y 2 F 1 (a + m, m; d; y) (A21) 

* — ' a, 77 77 L — ' c, 777 777 

n=l y 7 ' m=l v ' ; 



with 



which leads to 

T(b)F 4 (a, b; c, d; x, y)\ b=0 = T(0) + ^ 3 F a (l, 1, 1 + a; 2, 1 + d; y) + M {a, 0; c, d; x, y) (A22) 

a 

(1-y) 2 x 

Mq(<x, 0; c, d; x, y) = Go(a, 0, 1 — d + a, 1 — d; 1 — d + a, c; x, 1 — y) + L^la, b; c, d; , — ) — 

4x 4 

- r(d) y ~ ^ (1 - V) d - a 2F 1 (d,d-a;l-a + d;l-y). (A23) 

r(o) 

This result is not valid for y — 1, in which case we have to use 

r(6)F 4 (a, b; c, d; x, 1) | h=0 = T(0) + {ip{d) -ip(d-a))+ M (o, 0; c, d; x, 1) (A24) 
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with 

M (a, 0; c, d; x, 1) = G (a, 0, 1 - d + a, 1 - d; 1 - d+ a, c; x, 0) + Lo(a, 0; c, d; 0, -) (A25) 
and ^(x) the digamma function defined as r'(a;)/r(x). Additional auxiliary expressions are 

^ x m y n (a, m + n)T(m + n)(c, m)(d, m) 
G (a,0,c,d;e,/;a;,y) = > — — — ■ — — — r (A26) 



and 



Lo(o, 6; c, d; f ) = r(6)L(a, 6; c, d; f )• (A27) 



The structure of the regions of convergence is in this 
case more complicated and a more sophisticated algo- 
rithm to determine which one to choose is necessary. 



APPENDIX B: TENSOR BASES FOR THE 
THREE-GLUON VERTEX 

The most general basis for the three-gluon vertex con- 
sists of 14 tensors. They are build from all possible com- 
binations of independent external momenta and the met- 
ric tensor. Under certain circumstances this number re- 
duces. We list the used tensor bases here for reference. 

The simplest case is when the three external momenta 
fulfill 



2 2 3 2 

Pi = P2 = P3 = P > 



which leads to 



Pi ■ P2 = P2 ■ P3 = P3 ■ Pi = 



(Bl) 



(B2) 



This kinematic configuration is called the symmetric 
point and according to Celmaster and Gonsalves d9j only 
three tensors are necessary: 

T { j!l p {pi,P2,P3) = 

= H i(p 2 ) (Oi -P2) P 8 I1V + (p 2 -pa)nS vp + (p 3 -pi)J w ) - 



H 2 (P 2 ) 



(P2 ~P3)p{P3 -Pl)v(pi — Pa)p 



P 



tt I 2nP1 p ^2 m P3„ ~PUP2 P P3, 

H 3 {P ) 2 



(B3) 



To get the scalars we contract the integral with the basis 
tensors and get the matrix equation 



i=l 



Hi(P 2 ) T$i(pi,P2,p 3 ) rm Pl ,p 2 ,p 3 ), (B4) 



which can be solved for Hi, H 2 and H 3 . There are two 
ways to get numerical values for the contracted integrals: 



Either one decomposes all occurring scalars into p 2 , (q + 
Pi) 2 \ (<? ~ P2) 2 and q 2 and calculates the integrals, or 
one uses the method described below to calculate the 
complete vertex and contracts the result with the basis 
tensors. For the symmetric point we cross-checked both 
methods. 

The choice of the basis we used for the general kine- 
matic case was motivated by two things: First it simpli- 
fied calculations enormously, because it has its origin in 
the method described in |32j to reduce tensor integrals to 
scalar ones and needs only one or two integrals for each 
tensor instead of several. Second, because some of the 14 
original tensors coincide or the integrals are the same for 
different tensors, the number of necessary basis tensors 
decreases to 10: 

10 

T f?u P (Pi,P2,P 3 ) = ^2E. l (pi,p 2 ,p 3 ;K;d)T^ p (p 1 ,p 2 ,P3)- 

i=l 

(B5) 

These are 

r>i 7)1 r>i 

(B6a) 



1 , . Pl P PuPl P 

r^ P {Pi,P2) = -2 

Pi 

P2 p P2 u P2 p 
Pi 



pvp 



(Pl,P2) = 



(B6b) 



tL p (puP2) 



Pi 



(B6c) 



pvp 



(P11P2) 



Pl u P2 v P2 p +P2 u Pl u Pl +P2 u P2 vP l p 



Pi 



Tp]vp(puP2) = 9puPl p + 9ppPl v + 9vpPl^ 
{Pl,P2) = 9puP2 p + 9 W P2„ + gvpP2 fi , 



pvp 



rL p (pi,P2) 



pvp 



(PUP2) 



Pl P Pl,P2 p 

Pi 

Pl p P2 v P2 p 
Pi 



(B6d) 
(B6e) 
(B6f) 

(B6g) 



(B6h) 



Pl„P2 1 ,(P2 ~Pl)p + (P2 ~Pl)pPl v P2 p 



Pi -P2 



(B6i) 
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(B6j) 



We partially normalized them so that all of them have 
the same mass dimension. 

As this basis is the result of a mathematical method 
and not of physical considerations it does not have ex- 
plicit symmetries in the momenta. However, as we are 
only interested in direct numerical results we stick to this 
basis and do not transform to another one. The general 
formula for a tensor integral of rank M in ref. [32], eq. 
(11), has to be slightly changed, because it was originally 
derived for Minkowski space and we use the Euclidean 
metric: 



J ^ ) . MM ( fc i' • • • , fejv;^i, • •• , v N ; d) 



A,(Ti , . . . ,CTjv 

lA 



( _ 1) M (2) 2M-3A^M-A X 



X { [ 5 ] A [fc!p . . . [k N ]°» K *l) • • • (?N, °n) X 

x I {N) (pi, . . . ,p N ; vi + o-i, . . . , v N + a N ; d + 2(M - A)). 

(B7) 

This equation can be used for all tensor integrals of rank 
M with N legs as long as a formula for the scalar integral 
for arbitrary dimension d and arbitrary exponents z/j of 
the squared momenta appearing in the loop is at hand. 
For its definition in the case of the three-point integral 
see Eq. (|8}. For the use of Eq. (B7l a fixed momentum 
routing is necessary depicted in Fig. [14] To translate 
it to the momentum routing we have in Fig. [2] for the 
ghost triangle we set k% — 0, k 2 = ~P2 and k\ = p\. 
The ghost-ghost-gluon triangle has a different momen- 
tum routing and we set k 3 = and replace k 2 by p g h. 
The term in brackets is the symmetric tensor combina- 
tion of the metric tensor g and the momenta p\, . . ., p^. 
E.g. the tensor T3 stems from {p\p2} pup- The value of 
A goes from to the integer part of M/2 and the value 
of the <7j from to M. Thereby the sum is restricted 
such that only combinations of the metric tensor and the 
momenta occur which allow to distribute all M indices 
among them. The 16 terms we initially get for the ghost 
triangle can be combined to the ten tensors of Eq. (B6 1. 



We amended the tensors by adding normalization factors 
to get the correct canonical dimension and facilitate the 
comparison with the results of [H] . The technical advan- 
tage of this basis is the low number of integrals, namely 
14, necessary to be calculated. With the method de- 
scribed above for the symmetric configuration the num- 
ber of integrals necessary at least quadruplicates (using 
the basis where every basis tensor consists of only one 
expression), because the contraction of a tensor with the 
integrals yields four terms. Changing the basis increases 
the number of integrals due to more complex basis ten- 
sors. E.g. the Ball-Chiu basis 05] needs 46 x 4 = 184 
integrals, of which a few may have the same values for 
the exponents and thus coincide. The number of inte- 
grals is not important in the most part of the Euclidean 
regime, but at the boundaries convergence of the sums 
becomes slow so that less integrals mean quite a speed-up 
in calculations. 



k r k N 




k 5 -k 4 



FIG. 14: Initial momentum routing when using the method 
of ref. [32] for calculating tensor integrals. The Vi denote the 
powers of the squared momenta. 
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